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This  report  summarizes  the  work  done  under  Air  Force  Aero 
Propulsion  Laboratory  Contract  No.  F33615-78-C-2016  during  the 
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noise/performance  trade  procedure,  and  (2)  the  analytical  work 
describing  the  propagation  of  sound  waves  in  variable  area  soft 
walled  ducts  with  flow.  Ducts  of  this  type  are  found  in  the 
inlet  and  exhaust  passages  of  high-bypass  turbojet  engines. 

This  analytical  work  was  used  in  developing  computer  programs 
predicting  the  attenuation  of  sound  in  these  ducts.  The 
programs  are  detailed  in  user  manuals.  The  work  was  done  by  Dr. 
John  J.  Schauer,  Mr.  John  T.  Datko,  Mr.  Robert  W.  Guyton,  Mr. 
David  W.  Ewick,  Mr.  Eugene  P.  Hoffman,  and  Mr.  Robert  N.  Bergen 
of  the  University  of  Dayton  under  Project  Engineers  Mr.  Robert 
M.  McGregor  and  Dr.  James  S.  Petty  of  the  Aero  Propulsion 
Laboratory,  AFWAL/POTC,  Wright -Patterson  Air  Force  Base,  Ohio 
45433. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 

This  work  was  a  continuation  of  an  investigation  of  air¬ 
craft  duct  acoustics  which  started  at  the  University  of  Dayton 
in  1972.  In  this  period  the  University  of  Dayton  has  aided  the 
Air  Force  Aero  Propulsion  Laboratory  in  the  highly  analytical 
area  of  Duct  Acoustics. 

Under  this  most  recent  contract  the  effort  was  divided 
into  two  tasks.  Task  1  consisted  of  a  review  and  evaluation  of 
duct  acoustic  models  developed  by  others  and  published  prior  to 
or  during  this  effort.  Included  in  this  task  was  maintenance  of 
exiting  Air  Force  noise  programs  and  improvement  or  modification 
of  these  programs  in  the  area  of  duct  acoustics. 

Task  2  was  divided  into  three  phases.  Phase  I  consisted 
of  sound  transmission  prediction  in  ducts.  Of  particular 
concern  under  this  phase  were  the  examination  of  sheared  flow 
orthogonality  theories  previously  studied  and  the  variable  area 
duct  sound  transmission  problem.  Phase  II  required  the  inte¬ 
gration  of  previously  developed  models  into  the  primary  Air 
Force  noise  prediction  program.  Phase  III  ensured  that  the 
integrated  changes  worked  and  verified  the  resultant  program  by 
predicting  noise  for  various  duct  lining  configurations . 


TASK  1  -  RESULTS 

A  published  lining  design  model  was  adopted  as  a  part 
of  the  Propulsion  Lab  noise/performance  trade  methodology. 

The  resultant  program,  based  on  the  Boeing  publication  by 
Armstrong  tl],  is  an  improvement  over  the  previous  duct  lining 
program  because  it  includes  in  the  model  the  effects  of  boundary 
layer  thickness  and  local  sound  pressure  level.  It  also  models 
polymide  sheet  linings. 

Another  published  technique  was  developed  into  an  approxi¬ 
mate  method  for  estimating  sound  attenuation  in  lined  ducts. 

The  ray  acoustics  model  published  by  Rice,  Heidman  and  Sefrin 
1 2 ]  was  the  basis  for  this  program.  It  was  desirable  because  it 
produces  an  estimated  value  for  the  complex  axial  sound  propa¬ 
gation  constants  associated  with  the  duct  modes.  These  are 
required  as  an  initial  guess  by  the  eigenvalue  finding  program 
01 TSHE  developed  under  a  previous  contract.  Prior  to  this 
programs  development,  experience  was  the  only  basis  for  the 
initial  guess  required  for  OPTSHE  and  it  was  difficult  to  home 
in  on  the  desired  eigenvalue.  The  approximating  program  is  also 
useful  for  finding  a  quick  estimate  for  the  more  highly 
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attenuated  modes  which  play  only  a  minor  role  when  the  sound 
energy  is  distributed  among  all  the  propagating  modes.  This 
feature  was  utilized  in  the  Task  2  equal  energy  work. 


TASK  2  -  RESULTS 

As  part  of  Phase  I  of  this  task,  sound  transmission  in 
ducts  was  investigated.  In  particular  orthogonality  conditions 
in  shear  flow  and  transmission  in  variable  area  ducts  was 
studied. 

The  sheared  flow  orthogonality  proved  to  be  difficult. 
Interactions  of  the  acoustic  field  with  the  hydrodynamic  field 
were  evident.  Progress  in  the  uniform  mean  flow  situation  was 
good.  Here,  a  clearly  defined  fluctuating  hydrodynamic  field  was 
shown  to  fit  the  governing  equations  independent  of  the 
acoustical  field.  This  hydrodynamic  field  travels  down  the  duct 
at  the  speed  of  the  uniform  flow  velocity.  For  rotational  or 
sheared  flow  in  a  duct,  several  techniques  were  tried  to  produce 
a  solution  of  the  governing  partial  differential  equations.  No 
satisfactory  results  were  obtained. 

Variable  area  ducts  were  approached  with  a  transient 
finite  difference  solution  to  the  governing  equations.  The 
variable  duct  profile  was  mapped  conformally  into  a  rectangle. 
The  resultant  mapping  was  used  to  transform  the  governing 
equations  in  order  that  the  finite  difference  grid  could  be  set 
up  on  the  rectangle.  The  resuts  were  very  encouraging  for  the 
no  flow  conditions.  With  a  irrotation  mean  flow,  results  were 
also  obtained  but  not  over  the  range  of  parameters  of  interest. 

As  part  of  Phase  II  of  this  task,  the  theoretical  develop¬ 
ment  of  previous  contractual  efforts  were  integrated  into  the 
noise/performance  program.  Specifically,  the  least  attenuation 
mode  theory  and  the  equal  energy  theoretical  duct  attenuation 
were  included  as  options  to  the  existing  empirical  duct  attenua¬ 
tion  prediction.  In  addition,  the  program  was  expanded  to  permit 
the  results  of  any  attenuation  prediction  to  be  input  to  the 
program. 

Under  Phase  III  the  noise  prediction  program  output  was 
compared  to  flight  data  and  corrections  to  the  predicted  noise 
were  formulated  for  the  KC-135B  aircraft.  Hard  walled  engine 
ducts  and  three  soft  walled  ducts  were  evaluated  for  two 
different  engines,  the  TF33-P5  and  the  CFM  56.  The  physical 
properties  of  the  soft  walls  were  found  to  enable  the  lining's 
effect  on  performance  to  be  judged. 
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SECTION  II 


LINING  IMPEDANCE  MODELS 


PERFORATED  SHEET  IMPEDANCE  MODEL 

The  normalized  resistance  (R/pc)  and  normalized  reactance 
(X/pc)  can  be  calculated  for  a  perforated  sheet  acoustic  lining 
material,  using  Fortran  IV  Subroutines  ENGNOR  and  METNOR.  The 
former  works  entirely  in  english  units,  the  latter  exculsively 
in  metric  units.  (See  page  4  for  definitions  of  variables.) 

The  equations  used  for  the  calculation  of  both  the  norma¬ 
lized  resistance  and  reactance  were  adapted  from  Armstrong’s 
work.[l:p.  21]  Originally,  these  equations  combined  both 
english  and  metric  units,  however  they  were  modified  to  work  in 
uniform  units. 

The  following  are  the  equations  used  for  normalized 
resistance  in  both  ENGNOR  and  METNOR: 

r  ■ 

R/pc  =  RVDC  +  ZJ  +  -~==z  exp(SN)  (1) 


where  in  the  english  version  (ENGNOR): 

2 

I  *1*  f  K  ■  V<  I 

.  077t 


RVDC  = 


(Tt/519) 


s 


14.6 


OA 


(T. /519  +  .416) 


.0000383 (Tt/519)  *75/f  (t/d  +  1-OA) 


ZJ  = 


0A/ 14.7  (Tt 


/519  +  .416) 


(2) 


(3) 


E 


1.0251 


exp (- .5072t/d) 


(4) 


SN  =  -1.8 (2.54df (OA) /V) 2 


(5) 


The  RMS  total  particle  velocity  is  obtained  using: 


3 


DEFINITIONS  OF  VARIABLES  AND  UNITS 
IN  PERFORATED  SHEET  MODEL 

Units 


lymbol 

Definition 

English 

Metric 

c 

speed  of  sound 

ft/sec 

m/sec 

d 

perforate  hole  diameter 

in 

-1 

cm 

-1 

f 

frequency 

set 

sec 

L 

panel  core  depth 

in 

cm 

M 

mach  number 

— 

— 

OA 

percent  perforate  open 
area 

decimal 

percent 

P 

S 

static  pressure 

PSI 

N/m2 

t 

thickness  of  perforated 
sheet 

in 

cm 

Tt 

total  air  temperature 

°  R 

0  K 

V 

RMS  total  particle  velocity 

ft/sec 

cm/sec 

VSPL 

RMS  spectrum  particle 
velocity 

ft /sec 

cm/sec 

VGF 

RMS  grazing  flow  particle 
velocity 

ft /sec 

cm /sec 

0 

boundary  layer  momentum 
thickness 

in 

cm 

0 

local  air  density 

lbm/f t3 

kg/m3 

4 


~  *v 


my  i  iw jmw 


V  =  /v 


2  2 
SPL  VGF 


and  the  spectrum  particle  velocity  is: 


„  .  W144V 

SPL  pc 


where 


(  SPL 

PRMS  -  2’901  *  l0"9  \10  2°, 


The  grazing  flow  contribution  to  velocity  is: 


GF 


,25 (XKK)cM' 


RVDC  .  /RVDC2  .  XKK  „..2 

2~  V  2  +  —  EM 


where 


if  e  =  o 


XKK  = 


.05  +  . 11 (d/0)  if  0  /  0 
And  for  the  metric  version  (METNOR): 

2 


RVDC  = 


.  0  3031 5t 


101325 


OA 


(Tt/288) 
(Tt/288  +  .416) 


ZJ  = 


.0000383 (Tfc/288) ,72/f (T/d  +  1 


-  OA) 


“VTotIzS  <V288  +  -416) 


The  equation  for  E  is  the  same  as  equation  ( 4) 


SN  =  -1.8 (df (OA) /v) 


*.«.  -  ’-W  « 
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where 

w  -  W 

V__  *  (Eqn.  8)  x  100 

GF 


(12  ) 


The  normalized  reactance  was  calculated  in  ENGNOR  using: 

—  =  •Q0°46gf,r  {t  +  .85d(l-.7^OA)-{.0003657V/OA)  }  -  cot  (kL)  (13) 

PC  OA  y'Tt/519 

where 

k  =  2iTf/12c 
and  in  METNOR: 

x_  a  •-Q.O-QlMiL-  ft-.  4-  .aqdfl-.7/OA)-(.00003048V/QA)  }  -  cot(kL)  (14) 
Pc  OA  J  Tt/188 

where 

k  =  27rf/100c 


Subroutines  ENGHOL  and  METHOL  were  developed  to  solve  the 
previous  equations  for  the  percent  perforate  open  area  (OA)  and 
the  panel  core  depth  (L),  when  the  normalized  resistance  and 
reactance  are  given  quantities.  ENGHOL  is  the  english  version, 
while  METHOL  is  the  metric  version. 

Due  to  the  nature  of  the  equations,  the  solutions  for  (OA) 
utilizes  an  iterative  type  algorithm.  ENGHOL  calls  an  internal 
subroutine  LINPAR,  obtained  from  Schauer,  et.  al.[3:p.  142] 

( LINPAT  for  METHOL),  which  determines  an  initial  starting 
quantity  for  (OA),  which  is  then  used  in  the  iterative  scheme. 
The  appropriate  calculations  are  carried  out,  (R/pc)  is  calcu¬ 
lated  and  compared  to  the  input  value,  and  if  the  difference  is 
greater  than  .001,  a  new  value  is  calculated.  Once  (OA)  is 
determined,  (L)  is  solved  explicitly  from  equations  (13)  or 
(14),  depending  on  the  units  used. 

Polymide  Sheet  Impedance  Model 

The  normalized  resistance  (R/pc)  and  normalized  reactance 
(X/pc)  can  be  calculated  for  a  polymide  sheet  acoustic  lining 
material,  using  Fortran  IV  Subroutines  PIMENG  and  PIMMET. 

PIMENG  uses  only  english  units,  while  PIMMET  is  the  metric 
version . 
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DEFINITIONS  OF  VARIABLES  AND  UNITS 
IN  POLYMIDE  SHEET  MODEL 


Units 


Symbol 

Definition 

Enqlish 

Metric 

c 

speed  of  sound 

ft/sec 

m/sec 

f 

frequency 

-1 

sec 

sec  3 

L 

panel  core  depth 

in 

cm 

M 

mach  number 

— 

— 

Ts 

static  air  temperature 

°  R 

°  K 

V 

RMS  total  particle  velocity 

ft /sec 

cm/ sec 

VSPL 

RMS  spectrum  particle 
velocity 

ft/sec 

cm/sec 

VGF 

RMS  grazing  flow  particle 
velocity 

ft/sec 

cm/sec 

P 

local  air  density 

lbm/ft 3 

kg/m3 

N 

number  of  plies  (integer) 

— 

— 
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The  equations  used  in  the  calculations  were  obtained  and 
adapted  from  Armstrong . C 1 :p .  22] 

The  following  are  the  equations  used  to  obtain  the 
normalized  resistance  in  both  PIMENG  and  PIMMET. 

R/pc  =  (RVDC  +  S  +  T)/pC  ( 


where  in  PIMENG 


RVDC  =  4  60 • A'N1 ’ 333 


(16  ) 


and 


A 

S 

T 


1.115  x  10 

.  305-A-N2  *f 
2.86pN2V 


1.5 


,  T  +216 
L  s 


(17) 

(18) 


The  RMS  total  particle  velocity  is  obtained  using  equation  ( 6  ) , 
and  the  RMS  spectrum  particle  velocity  is  obtained  from  equation 
(  7  ).  The  grazing  flow  contribution  is  obtained  using 


.  05cM 


'GF 


RVDC 

pc 


RVDC' 


PC 


2  2 
. 71n  M^ 


(19) 


And  in  PIMMET 

1  333 

RVDC  =  22650-A-N 


and 


A  =  1.496  x  10 


-5 


rT  1.5 

s 

T  +120 
.  S 


S  =  1.4891 • A*N2 • f 
T  =  2 . 86pN2 (V/100) 


(20) 

(  21) 
(  22) 


The  RMS  total  particle  velocity  is  obtained  using  equation  (  6  ) 
and  the  RMS  spectrum  particle  velocity  is  obtained  again  using 
equation  (12).  The  TMS  grazing  flow  contribution  is: 

VGF  =  (Eqn .  19)  x  100  (23) 
The  normalized  reactance  is  obtained  in  PIMENG  using: 
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X 

pc 


[(.549pN1*5f) 


(5.12  x  10-5pN3fV) ] /pc  -  cot(kL)  (24) 


where 

k  =  2irf/12  c 
and  in  PIMMET  using: 

—  =  [( .1673pN1,5f)  -  (5.12  X  10“5pN3fV) ] /pc  -  cot(kL)(25) 
pc 

where 

k  =  2irf /100c 


Subroutines  LINENG  and  LINMET  were  developed  to  solve  the 
previous  equations  for  the  number  of  plies  (N)  and  the  panel 
core  depth  (L),  when  the  normalized  resistance  and  reactance  are 
given  quantities.  LINENG  is  the  english  version  while  LINMET  is 
the  metric  version. 

Because  of  the  nature  of  equation  (15),  the  number  of 
plies  required  to  obtain  a  given  resistance  cannot  be  solved  for 
explicitly.  Therefore  (N)  is  initially  given  a  value  of  1.0, 
then  equation  (15)  is  solved.  The  answer  is  compared  to  the 
input  value,  and  if  the  difference  was  greater  than  .001,  the 
value  of  N  was  modified  accordingly.  Once  N  is  determined,  (L) 
is  solved  for  explicitly  in  LINENG  using  equation  (24)  or  in 
LINMET  using  equation  (25). 

In  summary,  for  perforated  plate  or  polymide  sheet  duct 
linings  the  relationships  between  the  physical  dimensions  and 
the  acoustical  properties  have  been  programmed.  The  program 
gives  either  the  acoustical  properties  when  the  physical 
dimensions  are  input  or  vice  versa.  The  programming  is  required 
because  of  complexity  of  the  relationships  and  because  an 
iterative  solution  is  necessary  when  going  from  the  acoustical 
properties  to  the  physical  dimensions.  The  details  of  the 
program  are  given  in  the  program  users  manual. 
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SECTION  III 


AXIAL  PROPAGATION  CONSTANT  ESTIMATION 

For  a  cylindrical  duct  with  no  flow  we  would  like  to 
estimate  the  axial  propagation  constant,  x  -  i a,  for  a  given 
angular  and  radial  mode  number  (m,n)  and  for  a  given  wall 
impedance,  R  +  ix.  The  technique  is  to  estimate  r,  the  phase 
velocity,  for  the  desired  mode  by  assuming  a  hard  wall.  Then  to 
estimate  o,  the  duct  attenuation,  based  on  ray  acoustics  as 
given  by  Rice,  et  al  [2]. 

From  the  governing  partial  differential  equations  a 
separation  of  variables  gives  a  solution  in  terms  of  Bessel 
functions  in  the  radial  direction.  The  phase  velocity  for  a 
hard  wall  is  related  simply  to  the  zero's  of  the  derivative  of 
the  Bessel  function.  These  zero's  are  real  for  a  hard  wall  and 
as  can  be  seen  from  Rice,  et  al  [2]  equation  7  for  no  flow. 


(26) 


where 


Tmn 


+•  H  a-  V 

is  the  phase  velocity  of  the  m  angular,  n  radial 
mode 


n  is  the  frequency  parameter,  n  =  fD/c,  the  frequency 
times  the  duct  diameter  over  the  speed  of  sound 


and 


a 


n 


1.1.  i,  v 

is  the  n  ‘  zero  of  J  Ma  )  =  0,  J  being  the  nrn order 
Bessel  function.  m  n  m 


This  relationship  also  permits  identification  of  cutoff 

modes,  the  mode  which  does  not  propagate  in  the  hard  duct,  as 

the  modes  having  an's  large  enough  to  make  the  phase  velocity  x 

imaginary.  The  radial  propagation  angle,  4>r,  defined  by  Rice  to 

be  the  angle  of  incidence  of  the  wave  front  on  the  wall  is  given 

by  Rice  and  Heidman  [4],  equation  19  evaluated  at  the  wall  for 

real  a  as: 
n 


cos  <J>  =  —  / a  ^  (27] 
r  xr  rx  /  n 

The  axial  propagation  angle,  defined  by  Rice,  et  al  to  be 
the  angle  the  wave  front  makes  with  the  centerline  can  be  eval¬ 
uated  by  equation  5  of  Reference  2  as: 


Then,  consideration  of  a  wave  incident  on  the  wall  and 
reflected  back  with  the  angle  of  incidence  equal  to  the  angle  of 
reflection  and  meeting  the  soft  wall  boundary  condition  gives 
for  the  ratio  of  axial  energy  flux  after  the  bounce  to  that 
before  as  equation  25  of  Reference  2: 

I  (R2  +  X2)  cos2  <J>  —  2R  cos  4>  +  1 

_2L_  =  - - 2  - -  (28) 

Ix  (R  +  x  )  cos  4>r  +  2R  cos  4>r  +  1 

For  N  bounces  the  attenuation  constant  o  may  be  estimated 

as : 

„  _  ION  .  ^X_ 

°  ~  27. 3n  log10  +  (29) 

xx 

The  number  of  bounces,  N,  is  not  clear  because  of  the  var¬ 
iation  of  or  with  r  in  the  duct,  but  matching  the  estimated  value 
of  a  with  the  exact  value,  for  some  typical  cases  shows  that: 

1  J  an/lrT] 

N  =  i  tan  4>  =  '/  =  (30) 

/ 1  -  (an/7rn)  2 

gives  reasonable  agreement.  Equation  29  then  permits  the 
attenuation  constant  a  to  be  estimated  in  a  simple  direct  way. 
The  exact  technique  must  use  an  iterative  pro  ess,  with  an 
initial  guess  required  for  t  and  a.  Use  of  these  estimates  as 
the  initial  guess  usually  permits  the  exact  technique  to 
converge  to  the  desired  mode. 

These  relationships  were  programmed  to  permit  their 
repeated  evaluation  for  all  propagating  modes.  The  details  of 
the  program  are  in  the  Users  Manual  for  Suppressor  Performance 
Estimation  included  as  part  of  this  report.  The  Users  Manual 
compares  the  estimates  to  the  exact  values  and  shows  that  the 
estimates  are  poor  when  the  radial  propagation  angle  is  close  to 
zero  or  V2  radians.  This  is  to  be  expected  since  these  angles 
correspond  to  no  bounce  or  infinitely  many  bounces  as  the  wave 
propagates  down  the  duct.  Infinitely  many  bounces  occur  as 
cutoff  is  approached.  Zero  bounces  corresponds  to  the  hard  wall 
propagation  of  a  plane  wave.  For  most  of  the  higher  order  modes 
which  are  not  too  close  to  cutoff,  however,  the  procedure  gives 
useful  results. 


SECTION  IV 


SHEARED  FLOW  ORTHOGONALITY 


INTRODUCTION 

An  orthogonal  condition  which  permitted  isolation  of  the 
modal  coefficients  for  an  expansion  of  an  arbitrary  initial 
pressure  distribution  in  terms  of  the  sheared  flow  pressure 
modes  was  investigated  in  an  Aero  Propulsion  Laboratory  report 
[3].  In  addition  to  a  satisfactory  solution  of  the  adjoint 
problem  for  the  second  order  partial  differential  equation,  it 
was  recognized  that  the  second  order  equation  did  not  permit  an 
arbitrary  initial  velocity  both  axially  and  radially,  as  is 
physically  possible.  This  led  to  the  effort  reported  here,  which 
includes  investigations  in  several  areas.  Hydrodynamic  modes 
were  investigated  in  uniform  flow  and  sheared  flow.  Because  of 
difficulty  with  the  sheared  flow  hydrodynamic  modes,  a  simple 
linear  velocity  profile  flow  with  a  small  Mach  number  was 
studied.  The  perturbation  approach  used  did  not  result  in  a 
solution  which  satisfied  all  the  boundary  conditions. 


HYDRODYNAMIC  MODES  IN  UNIFORM  FLOW 

For  a  rectangular  duct  with  uniform  flow  M  in  the  axial,  or 
x  direction  the  continuity,  x  and  y  momentum  equations  for  the 
perturbations  are  respectively: 


u  +  v  -  - 
x  y 


pc 


7  [pt  +  CMV 


u.  +  cMu  =  -  —  P 
t  x  p  X 


(31) 


vt  +  CMVX  -  -  -  Py 

where  u  and  v  are  the  axial  and  cross  duct  velocities  and  c  is 
the  speed  of  sound.  Here  p  is  the  mean  density.  For  sinusoidal 
motion  in  time  u,  v,  and  P  are  proportional  to  exp(iwt)  and  for 
k  =  u>/c  we  get : 

“x  +  vu  -  *  js  (lkp  +  “V 


iku  +  MUX  -  -  -  Px 


(32) 


ikv  +  Mvx  =  -  —  Py 
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Looking  for  a  hydrodynamic  solution  which  has  no  fluctuating 
pressure  associated  with  it  would  mean  setting  P,  px ,  p  equal 
to  zero  in  Equation  32  leaving  y 

u  +  v  =  0,  iku  +  Mu  =  0,  ikv  +  Mu  =  0  (33) 

x  y  *  x 

multiplying  the  momentum  equations  (Equation  33)  by  dx  and 
integrating  while  holding  y  constant  gives: 

u  =  e  M  •  f  (y)  and  v  =  e  M  •  g(y) 

where  f  and  g  are  arbitrary  functions  of  y  except  that  the 
boundary  condition  on  v  at  the  duct  walls  is  that  v  must  be  zero 
in  order  to  avoid  effecting  the  acoustical  wall  impedance.  To 
meet  this  boundary  condition  on  v  at  y  =  0  and  at  y  =  1,  g(y)  can 
be  built  up  by 

=<  - 

v  =  e  J  n»y  (35) 


Substituting  Equation  (35)  into  the  continuity  equation  of 
Equation  (33)  gives: 


u 


00 


I 

n=l 


bnsin  mry 


(36) 


where  an  and  must  be  related  by: 


a 

n 


-iis  b 

nirM  n 


(37) 


Further,  taking  (  )x  of  the  y  momentum  equation  and  (  )y  of  the 
x  momentum  equation  of  Equation  (33),  subtracting  and  solving  the 
resultant  equation  as  before  gives  the  fluctuating  vorticity  as 


v  -  u  =  e 
x  y 


ik 

M 


h(y) 

where  h(y)  is  related  to  f  and  g  such  that: 

sin  mry 


(38) 


h(y)  =  l 

n=l  L 


k)2  ^  .  .2 

+  <n7r> 


1  — 
J  nit 


(39) 


The  fluctuating  vorticity  then  travels  down  the  duct  in 
cells  of  length  Mc/f ,  where  f  is  the  frequency  and  moving  at  the 
speed  of  the  flow.  Me,  as  can  be  seen  by  putting  the  time  depen¬ 
dence  back  in,  giving  u,  v,  and  vx  -  Uy  all  proportional  to 
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exp(-ik/M( x-Mct )) .  These  fluctuating  velocities  without  fluc¬ 
tuating  pressures  are  called  the  hydrodynamic  modes. 

A  similar  solution  can  be  found  in  axisymmetric  uniform 
mean  flow  situations. 


4 
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SHEARED  PLOW  STUDY 


For  a  rectangular  duct  with  a  sheared  flow  where  M  is  a 
function  of  y  only,  the  continuity,  x  and  y  momentum  equations 
with  sinusoidal  time  dependence  for  the  perturbations  are 
respectively: 


ux  +  v 
iku  + 

ikv  + 


y 

Mu 

Mv 


x 


X 


+ 


1_ 

DC 


(ikP  +  MPx) 


1_ 

pc 


P 

y 


(40) 


These  differ  from  the  uniform  flow  equations  (Equation  32)  only 
by  the  addition  of  a  single  term  in  the  x  momentum  equation. 
Although  this  seems  to  be  only  a  small  difference,  it  causes  a 
coupling  between  the  fluctuating  velocity  field  and  the  fluc¬ 
tuating  pressure  field.  While  the  irrotational  mean  flow 
equations  (Equation  32)  can  be  reduced  to  a  second  order  equation 
in  pressure  alone  these  sheared  flow  equations  reduce  to  a  third 
order  partial  differential  equation  in  pressure  alone. 

Although  considerable  effort  was  expended  in  an  effort  to 
find  solutions  to  the  set  of  equations  (Equation  40),  no 
solutions  matching  duct  boundary  conditions  were  found.  Among 
the  techniques  tried  were  perturbation-type  solutions  for 
M  =  Mq  +  ey,  where  e  is  a  small  parameter. 
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SECTION  V 


VARIABLE  AREA  DUCTS 


INTRODUCTION 

Because  aircraft  ducts  have  a  variable  cross-sectional 
area,  the  sound  propagation  in  such  ducts  is  of  interest.  The 
complexity  of  the  governing  partial  differential  equation  and 
boundary  conditions  for  fluctuating  pressures  and  velocities 
under  these  conditions  has  led  to  numerical  methods  as  the 
appropriate  technique  for  solving  these  governing  equations.  A 
November  1980  survey  of  numerical  techniques  by  Baumeister  [5] 
shows  the  progress  in  this  area  and  the  current  limitations. 

This  effort  is  in  good  agreement  with  Baumeister' s  comments  that 
transient  numerical  techniques  which  step  the  governing  equa¬ 
tions  ahead  in  time  are  very  promising.  The  present  effort 
combines  these  transient  techniques  as  detailed  by  Baumeister  in 
References  6  and  7  with  mapping  techniques  introduced  by  Quinn 
[8],  Quinn's  analytical  maps  for  particular  duct  shapes  are 
extended  here  to  a  numerical  mapping  technique  for  general  duct 
shapes  including  bullet  nose  type  aircraft  inlets.  The  mapping 
technique  consists  of  a  coordinate  transformation  which  maps  the 
duct  boundaries  into  a  actangle.  The  governing  partial  differ¬ 
ential  equations  are  transformed  by  this  coordinate  change  into 
a  new  set  of  partial  differential  equations  which  may  be  solved 
together  with  the  transformed  boundary  conditions  in  the  rec¬ 
tangle  using  finite-difference  methods. 

Two  programs  were  developed.  The  first  program  called 
TDEP  solved  the  no  flow  problem  in  a  variable  area  duct  for 
arbitrary  wall  impedance  distribution.  The  second  program, 
called  TDEPF ,  solved  the  governing  equations  with  flow  for  an 
axisymmetric  potential  flow  distribution.  This  irrotational 
flow  gives  no  coupling  between  the  pressure  distribution  and 
hydrodynamic  modes . 


NO  FLOW  VARIABLE  AREA  DUCT  ACOUSTICS 


The  asixymmetric  no  flow  governing  partial  differential 
equations  can  be  reduced  to: 


P 


z  z 


+ 


P  +  -  P  + 
rr  r  r 


(41) 


as  given  for  example  by  Baumeister  [7],  Here  r,  z,  and  9  are 
the  usual  cylindrical  coordinates  with  lengths  nondimensiona- 
lized  by  dividing  by  R,  P  is  the  fluctuating  pressure  non¬ 
dimens  ionali  zed  by  pc2  and  nr  =  Rf/c,  the  initial  duct  radius  or 


annular  radius  difference  times  the  frequency  f,  divided  by  the 
acoustic  velocity  c.  The  mean  density  p  is  used  in  making  P 
dimensionless.  Time  is  also  made  dimensionless  here  by 
t  =  f • t ime . 

The  9  variation  is  assumed  to  be  proportional  to  e^1110 
giving  Pee  =  -m2P.  Equation  (41)  then  is  reduced  to  two  space 
variables.  Following  Quinn  [83  we  define  a  coordinate  transfor¬ 
mation  z  =  f(x,y),  r  =  g(x,y).  Then  for  a  conformal  coordinate 
transformation  satisfying  the  Cauchy-Riemann  conditions  fx  =  gy 
and  fy  =  -gx  Equation  (16)  transforms  to: 


p  +  P  +  h  Py  '  h  P*  -  ^  <fx2+fy2)P  '  (fx2+fy2>''2ptt  1421 


xx  yy 


The  boundary  conditions  on  a  soft  wall  C  =  P/V',  where  C 
is  the  wall  impedance  over  pc  and  V'  is  the  normal  component  of 
velocity,  are  substituded  into  the  normal  momentum  equation: 


(43) 


to  give 


(44) 


as  the  boundary  condition  in  the  duct.  This  transforms  to: 


n (f  tf  2>1/2 

P  =  -  1_2 - £ -  P 

y  c,  t 


(45) 


for  the  outer  soft  wall  boundary  condition. 


In  a  similar  manner  the  exit  impedance  c,e  appears  in  the 
duct  exit  boundary  condition  as: 


P  =  - 
x 


n  (f  2+f  2)1/2 
1  x  y _  p 

^  t 


(46) 


The  exit  impedance  was  taken  from  Baumeister  [7]  to  be 

-  II  -  I'W'S’V2!172- 

The  inner  wall  of  annular  sections  was  taken  as  hard  but 
it  could  easily  satisfy  a  soft  boundary  condition  Equation  (45) 
with  the  sign  on  Py  changed  to  negative.  On  the  centerline  of 
cylindrical  ducts  the  conditions  applied  were  Py  =  0  for  angular 
mode  m  =  0  or  P  =  0  for  m  >  1. 
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The  entrance  condition  used  follows  Baumeister  [7]: 

P  (o,y,t)  =  (47) 

where  Jm  and  were  as  in  Equation  (26) 

The  initial  condition  in  the  duct  was 

P  (x,y ,o)  =  0;  Pt(x,y,o)  =  0  (48) 


The  finite  difference  representation  of  Equation  (42)  and 
its  boundary  conditions  was  based  on  central  differences,  again 
following  Baumeister  [7].  Table  1  shows  the  coefficients  in: 


where 


p^+l  =  2Pk  .  -  p*"1  +  - L _  ju 

i,j  1,3  1,3  n2(4y)2(£x2+fy2)  1 


h  -  vti.j +  +  vli.j 


+  E  Pk  +  P  pk'*"  ^  +  G  pk— ^ 

i  i+l(j  4  i,j  4  i,j 


(49) 


(50) 


and  the  £  subscript  denotes  the  cell  location.  Cell  1  is  an 
interior  cell,  cell  2  on  the  outer  boundary,  cell  3  on  the  inner 
boundary  or  centerline,  cell  4  in  the  exit  plane,  cell  5  the 
outer  boundary  exit  plane  corner  cell,  and  cell  6  the  inner 
boundary  exit  plane  corner  cell. 


As  seen  from  Table  1  the  coefficients  in  Equation  (50) 
depend  on  the  partial  derivatives  of  the  mapping  function 
2  =  f(x,y)  and  on  g  from  r  =  g(x,y).  These  were  obtained  from 
the  conformal  mapping  technique  as  detailed  under  the  subheading 
"Duct  Mapping". 


The  steady-state  pressures  were  evaluated  from  P(x,y)  = 
P(x,y,t)/exp(  2-n-it) .  Stability  of  the  system  was  greatly 
improved  by  adding  damping  with: 


,k+  1 


=  Pk+1  + 


DD  [Pk-1e4 


i  t 


-  Pk+1] 


(51) 


at  each  grid  point  after  a  pass  to  compute  new  pressures.  Here 
DD  is  the  damping  factor  with  no  effect  for  DD  =  0  and  no  change 
in  the  steady-state  solution  for  DD  =  1.  A  DD  of  0.2  seemed  to 
damp  oscillations  well  without  changing  the  steady-state  solu¬ 
tion.  The  addition  of  this  damping  seemed  to  permit  soft  walled 
solutions  to  be  stable  with  a  hard  walled  initial  pressure  mode 
at  the  duct  entrance. 

i 

1 
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Baumeister's  straight  duct  no-flow  results  of  Reference  7 
for  the  third  angular  mode  were  duplicated.  Quinn's  analytical 
conical  duct  solution  was  compared  to  the  TDEP  solution  for  the 
zero  angular  mode  giving  a  good  comparison  with  damping  intro¬ 
duced.  Quinn's  hard  wall  hyperbolic  horn  for  the  third  angular 
mode  was  also  run  as  a  test  case  for  the  acoustic  transient 
technique.  The  zero  angular  mode  in  a  soft  walled  straight 
duct  approached  the  same  steady  solution  with  several  different 
values  of  damping  introduced.  Several  variable  area  duct  con¬ 
figurations  were  run  satisfactorily.  The  bullet  nose  inlet 
configuration  is  shown  in  Figure  1.  The  resultant  pressure 
distribution  along  the  bullet  nose  duct  outer  wall  surface  for  a 
plane  wave  is  shown  in  Figure  2.  These  results  are  for  a  dimen¬ 
sionless  wall  impedance  5  =  9  -  i  on  the  outer  wall  with  the 
inlet  bullet  hard. 


INCOMPRESSIBLE  FLOW  VARIABLE  AREA  DUCT  ACOUSTICS 

Writing  all  the  variables  in  dimensionless  form  using  the 
initial  duct  radius,  R,  the  sound  frequency,  f,  the  speed  of 
sound,  c,  and  the  mean  density,  p,  gives  r  =  r*/R,  z  =  z*/R, 
t  =  ft*,  P  =  P*/  c  ,  u  =  u*/c,  v  =  v*/c,  w  =  w*/c,  U  =  U*/c, 

V  =  V*/c,  W  =  W*/c.  Here  the  stared  variables  refer  to  the 
dimensional  quantities.  With  n  =  Rf/c  the  following  set  of 
partial  differential  equations  governs  the  fluctuating 
variables : 
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These  equations  representing  continuity  and  the  three 
momentum  equations  are  analogous  to  those  given  by  Quinn  1 9 ]  and 
Eversman,  et  al  110].  In  the  last  equation  of  the  set  (Equation 
(52)  the  assumption  that  P  is  proportional  to  exp(ime)  has  again 
been  incorporated.  With  our  coordinate  transformation 
z  =  f(x,y),  r  =  g(x,y)  satisfying  the  Cauchy-Riemann  conditions. 
Equations  (52)  transform  to: 

f2nPt  -  f yVx  +  fxVy  +  £XUX  +  fyuy  +  f 2 |V  +  imw) /c< 


+  P(-f  V  + 
y  x 


f  V  + 
x  y 


f  u  + 
x  x 


f  u  + 

u  y 


f2V/g) 


+  (f  U  -  f  V)P  +  (f  U  +  f  V)P  =  0 
x  y  x  y  x  y 


19 


Normalized  Axial  Distance,  Z/R 


Figure  1 


J 


Bullet  Nose  Duct  Profile 


<U 

S-l 

3 

(0 

UJ 

<D 

a. 


(0 

O 

« 


Normalized  D] 


Figure  2.  Pressure  Distri 


2 


f2nut  +  (£XU  -  fyV)ux  +  tfyU  +  £xV)uy 
+  (£XU  -  fyv  +  fxPU  -  f  PV)0x 

+  (fyU  +  fxV  +  fyPU  +  f^PVlUy  +  £XPX  ♦  fyPy  =  0 


f2nvt  +  (fxu  *  fyv)v*  +  (£yU  +  fxv)vy 

+  (fu-fv+fPU-f  PV)V 
x  y  x  y  '  x 

+  (fyu  +  fxv  +  fyPU  +  fxPV)Vy  -  fyPx  +  fxPy  =  0 
f2nwt  =  (fxU  ~  fyV)wx  +  {fyU  +  fxV)wy  +  f 2 (w  +  imP) =  0 


where  f2  =  fx2  +  fy2. 

The  boundary  conditions  on  the  soft  outer  wall  of  the  duct 
should  be: 


V  =  Pt/c  +  U<  {P/c)  t/ti 


where  V  is  the  normal  velocity  component,  U'  is  the  mean 
tangential  velocity  component  and  x  denotes  the  tangential 
direction . 

Because  of  stability  problems  the  boundary  condition  used 
was  simply: 
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p  =  xy'  (55) 

with  u,  v  and  w  computed  from  the  momentum  equation,  the  bottom 
three  equations  in  Equation  (53) . 

On  the  inner  wall  or  centerline  the  boundary  conditions 
used  were: 


Pjj  =  0  and  v '  -  0 ,  w  =  0 
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where  Pn  is  the  normal  derivative  of  pressure  and  the  z  direc¬ 
tion  velocity,  u,  is  computed  from  the  z  momentum  equation,  the 
second  equation  in  Equation  (53) . 

The  inlet  boundary  condition  was  taken  as: 

P  =  J  (a  y)exp(2rrit)  (57) 

m  mn 

where  amn  was  the  radial  engenvalue  for  a  cylindrical  duct  with 
the  inlet  cross-section.  The  variables  u,  v  and  w  were  computed 
from  the  momentum  equations  of  Equation  (53)  at  the  inlet. 

The  exit  plane  was  computed  from  P  =  CeUj^  where  un  is  the 
normal  velocity  component  and  Ce  is  taken  from  Baumiester  [7] 


C  =  (1  -  M8)/(B-M) 


where 


B  =  Ul  -  (a  /27rn)2]  (1  -  M2)  > 
mn 


Again  u,  v  and  w  were  computed  from  the  momentum  equation  in 
Equation  (53) . 

The  initial  condition  in  the  duct  was  taken  as  P=u=v=w=0 
or,  in  some  computer  runs,  the  steady  solution  was  estimated  and 
used  as  an  initial  condition  to  reduce  run  time. 

The  finite  difference  representation  of  Equation  (53)  and 
the  boundary  conditions  was  based  on  a  predictor-corrector 
MacCormack  method  with  generalization  as  suggested  by  Petty  [11] 
to  cause  truncation  errors  in  both  to  be  proportional  to  the 
space  difference  squared  (second  order  accurate).  The  top 
equation  of  Equation  (53)  ,  the  fluctuating  continuity  equation 
was  written: 

k+T  k 

P.  .  =  P.  .-(At/nf.,)  (space  derivatives  in 

13  13  2  v.u.w.U.V,  and  P,  (6°> 


for  a  predictor  and 


=  (P*j  +  P^1)/2 


(same  space 
derivatives) 
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as  a  corrector.  The  space  derivatives  specified  by  Equation 
(53)  were  evaluated  in  the  predictor  using  a  forward  difference 
at  interior  points  and  where  it  was  feasible  on  the  boundary. 

At  boundary  points  where  forward  differences  were  not  feasible 
backward  differences  were  used.  In  the  corrector  the  space 
derivatives  were  evaluated  by  a  subroutine  which  used  backward 
differences  at  interior  points  and  three  point  formulations  as 
developed  by  Petty  [11]  to  maintain  the  second  order  accuracy  at 
the  boundary  points. 

The  time  and  space  derivatives  were  handled  in  an 
analogous  manner  for  the  time  derivatives  of  u,  v  and  w  in  the 
fluctuating  z,  r  and  8  momentum  equations  respectively  of  Equa¬ 
tion  (53)  . 

Where  boundary  conditions  were  specified  in  Equations  (54) 
through  (58)  the  appropriate  equation  from  Equation  (53)  was 
replaced  by  the  boundary  condition  equation  and  the  other  three 
quantities  were  stepped  forward  in  time  using  the  remaining 
three  equations  of  Equation  (53)  . 

Some  difficulty  was  experienced  in  obtaining  stable 
solutions  of  this  set  of  equations  and  boundary  conditions. 
Damping  the  solution  by  equations  of  the  form  in  Equation  (51) 
did  not  help  stablize  the  solution.  Different  combinations  of 
boundary  conditions  were  tried  but  no  formulation  was  found 
which  produced  stable  results  over  the  entire  range  of  para¬ 
meters  desired.  The  solution  seemed  sensitive  to  the  specifi¬ 
cation  of  the  boundary  conditions.  The  formulation  specified 
produced  good  solutions  for  a  slightly  flared  duct  at  Mach 
numbers  of  -.2  to  +.1  as  shown  in  Figure  3. 


DUCT  MAPPING 

The  variable  area  duct  acoustics  problem  was  solved  by 
mapping  the  interior  of  the  duct  in  terms  of  z  and  r  into  the 
interior  of  a  rectangular  region  specified  by  x  and  y.  The 
mapping  used  was  based  on  the  conformal  Schwartz-Chr istof fel 
transformation  as  presented  in  Reddick  and  Miller  [12]  (page 
423).  Following  the  standard  technique  the  duct  shown  in  Figure 
1  was  mapped  via: 


dw  _  „,_ttF  _  x  ai /u-1  ,^irF _ x^/tt-I 


WIT*  */  /  _  '  i  \ 

dF  =  K(e  _al) 


(e  -a2) 


(e7rF-a6)  as /lT_1 


(62) 


Here  w  *  z  +  ir  and  F  *  x  +  iy.  The  real  constants  ai  are 
interpreted  as  vertex  locations  on  the  real  axis  in  an  inter¬ 
mediate  plane,  the  Z  plane  where  Z  *  exp(itF).  The  duct  was 


Outer  Mall 
Imaainarv  Pressure 


Outer  Wall  Outer  Wall 

Imaginary  Pressure  Real  Pressure 


assumed  to  have  zero  included  angle  between  the  sides  at  the 
inlet  and  differential  in  radius  of  1  at  the  inlet.  That  is.  Aw 
at  Z  =  0  or  F  of  minus  infinity  is  i.  Expressing  w  in  expo¬ 
nential  form  and  comparing  the  change  in  w  equal  i  for  Z  going 
from  R  exp(O)  to  R  exp(iri)  as  R  approaches  zero  gives  K  as: 

K  -  l/7r-(-a1)ai/ff':L-{-a2)a2/lT‘1*--(-a6)a6/'TT'1  (63) 


The  negative  values  of  a^  corresponding  to  the  outside  wall  are 
straightforward  in  Equation  (63)  and  the  positive  values  corres¬ 
ponding  to  the  inner  wall  which  result  in  a  negative  value  to  a 
power  are  evaluated  from: 

(-ai)  1  =  [  1  a±  1  *exp(iri)  ]  1  =  (aj  ^-expfTri^)  (64) 

In  order  to  use  the  mapping,  the  radius  at  a  given  value  of 
x  and  y  must  be  found  because  it  appears  directly  in  the  govern¬ 
ing  partial  differential  equations.  The  values  of  w  =  z  +  ir 
are  found  from  the  integration: 

w  =  J  ||  dF  =  /  ||  d(x  +  iy)  (65) 


or 


“  -  /  3?  d*  +  *  f? 


(66) 


form: 


For  the  dw/dF  used  (Equation  62)  this  integral  is  of  the 


w  =  J  Pdx  +  Qdy 


(67) 


with  3P/3y  =  3Q/9x  except  on  y  =  0.  Reddick  and  Miller  [121 
(page  441)  shows  that  Green's  lemma  implies  that  the  integral  of 
Equation  (65)  is  independent  of  path  for  these  conditions  as 
long  as  all  paths  are  above  the  singularities  on  y  =*  0.  A 
numerical  integration  in  the  F  plane  was  used  to  evaluate 
w  »  z  +  ir  at  all  grid  points.  The  integration  was  carried  out 
by  integrating  horizontally  along  the  j  *  4  set  of  grid  points 
and  then  integrating  vertically  for  the  various  i  values. 
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With  the  conformal  mapping  established  it  can  be  shown 

that : 

'  >*  -  'V  >x  +  V  'y]'(£x2+  fy2)  <«8> 

and 

<  >r  -  'V  >y  -  V  »xJ/<fx2+  fy2) 

These  relationships  permit  transformation  of  the  governing 
partial  differential  equations  (Equation  41  or  52)  and  their 
corresponding  boundary  conditions. 


SECTION  VI 

NOISE/PERFORMANCE  MODEL  INTEGRATION 


LEAST  ATTENUATED  MODE  INTEGRATION 

Previous  contractual  efforts  had  produced  a  theory  for  the 
maximum  attenuation  of  the  mode  attenuated  least  in  rectangular, 
annular,  and  cylindrical  ducts.  This  theory  was  for  constant 
area  ducts  with  flows  at  Mach  numbers  of  magnitude  less  than  0.5 
and  boundary  layer  thicknesses  less  than  seven  percent  of  the 
duct  radius.  The  integration  of  this  work  into  the  noise  model 
program  developed  by  Dunn  and  Peart  [13]  was  simplified  by  the 
conclusion  that  the  maximum  attentuation  was  independent  of  the 
boundary  layer  thickness  over  this  range.  The  optimum  wall 
lining  as  shown  by  the  previous  efforts  in  Reference  14  was 
sensitive  to  the  boundary  layer  thickness. 

A  study  of  the  numerical  results  of  this  previous  effort 
permitted  relationships  between  peak  attenuation  and  frequency 
parameter  to  be  established  for  cylindrical,  annular,  and  rec¬ 
tangular  ducts.  These  relationships  were  put  in  arrays  as  an 
alternative  to  the  empirical  arrays  in  the  original  program.  As 
shown  for  the  cylindrical  duct  in  Reference  3,  page  39,  the 
least  attenuated  mode  theory  predicted  a  larger  attenuation  at 
low  frequency  than  Dunn  and  Peart's  empirical  predictions  and 
predicted  a  lower  attentuation  at  high  frequencies  than  the 
empirical.  The  difference  at  low  frequencies  was  possibly 
associated  with  deviations  from  the  assumed  conditions  which 
would  cause  measured  attenuation  to  be  less  than  that  predicted. 
At  high  frequencies  the  attenuation  predicted  by  the  least 
attenuated  mode  theory  was  lower  than  the  empirical  because  the 
theory  was  very  conservative  when  a  large  number  of  modes 
propagated  in  the  duct  as  happens  at  high  frequencies.  All  of 
the  other  modes  were  attenuated  more  than  the  least  attenuated 
mode  on  which  the  prediction  was  based. 

For  these  reasons  the  least  attenuated  mode  theory  was 
suspect  at  high  frequencies  and  other  assumptions  regarding  how 
much  of  the  sound  energy  was  in  each  of  the  propagating  modes 
were  suggested. 


EQUAL  ENERGY  THEORY  INTEGRATION 

The  equal  energy  theory  assumed  that  an  equal  portion  of 
the  total  sound  energy  was  in  each  propagating  mode.  Propagat¬ 
ing  modes  were  those  with  real  phase  velocities  as  defined  in 
Equation  (26)  .  Since  the  a  in  Equation  (26)  was  for  hard  walls 
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there  was  the  possibility  that  another  mode  or  two  might  propa¬ 
gate  at  high  frequencies  but  the  modes  which  barely  propagated 
were  highly  attenuated  and  did  not  effect  the  results.  At  the 
highest  value  of  the  frequency  parameter  in  the  original  program 
range  there  were  114  modes  propagating. 

Rather  than  finding  the  exact  attenuation  of  every  mode 
with  the  usual  iterative  process,  the  results  of  Section  III 
were  used  to  estimate  the  modal  attenuations.  These  estimates 
were  used  to  start  the  iterative  process  for  the  lightly 
attenuated  modes  which  contribute  most  to  final  energy  at  the 
end  of  the  duct.  The  attenuations  of  more  heavily  attenuated 
modes  which  have  only  a  small  contribution  to  the  final  energy 
were  taken  directly  from  the  approximate  results  of  the  Section 
III  analysis. 

The  attenuation  of  these  modes  with  equal  energy  in  each 
at  the  duct  inlet  was  calculated  from 

.  N  -kLo . 

VEi  -  5,1,'  <69> 


or  in  decibels. 

Attenuation  =  10  log.,.  E  /E .  (70) 

310  o  i 

Here  E0  and  Ei  are  the  energies  out  and  in  the  duct.  The  number 
of  propagating  modes  is  N  and  the  attentuation  constant  is 
found  from  Equation  (29)  or  an  iterative  solution.  The  reduced 
frequency  based  on  duct  length  kL  also  appears. 

From  these  results  the  predicted  peak  attenuation  as  a 
function  of  frequency  parameter  and  normalized  duct  length,  L/H, 
was  found.  An  array  containing  this  information  was  included  in 
the  noise  prediction  program  as  a  second  alternative  to  the 
empirical  prediction  of  peak  attenuation.  This  array  of  peak 
attenuation  values  as  a  function  of  y  =  HF/c  and  L/H  is  shown  in 
Table  2.  This  table  is  used  only  under  control  option  IMAm  =  3 
as  detailed  in  the  subsection  entitled,  "Modification  Control". 


EXTERNAL  ATTENUATION  MODIFICATION 

The  preceding  two  theories  are  by  no  means  the  only 
possible  theoretical  attenuation  approaches.  They  do  not 
consider  variations  of  the  duct  lining  along  its  length.  They 
do  not  consider  variations  of  duct  cross-sectional  area.  They 
do  not  consider  arbitrary  distribution  of  the  sound  energy  among 
the  propagating  modes.  Numerical  schemes  such  as  the  finite 
difference  approaches  in  Section  V  can  include  these  effects. 
Consequently,  the  noise  prediction  program  was  modified  to 
permit  the  results  of  an  external  program  like  those  of  Section 
V  to  be  used  in  the  noise  prediction. 
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The  approach  to  this  modification  was  simply  to  bypass  the 
internal  attenuation  process  and  to  enter  the  externally  com¬ 
puted  attenuation  of  the  desired  duct  directly  as  attenuation  as 
a  function  of  frequency.  This  was  accomplished  using  existing 
control  parameters  and  simply  changing  their  range. 


MODIFICATION  CONTROL 

The  existing  parameters  LGMm  and  IMAm  were  used  to  control 
the  new  program  modifications.  As  they  existed  prior  to 
modification  in  Crawley,  et  al  C15]  the  old  usage  was: 


OLD  PARAMETER  USAGE 

IMAm  Specifies  whether  program  calculates  or  used  defines  peak 
attenuation  for  each  target  frequency. 

=  0  program  calculates 

=  1  user  defines  PLAm  values 

LGMm  Specifies  whether  program  calculates  peak  attenuation 
using  lining  geometry  or  user.  Defined  effective  duct 
height  and  ratio  of  treatment  length  to  effective  duct 
height . 

*  0  user  inputs  EDHm  and  ELOHm 

=  1  user  inputs  RlWm  and  TLm 

(the  =  1  option  only  if  IMAm  =  0) 

Here  PLAm,  EDHm,  ELOm,  RlWm,  and  TLm  are  defined  in  Reference 
[15].  As  currently  modified  the  new  usage  of  these  parameters 
is : 


NEW  PARAMETER  USAGE 

IMAm  Specifies  how  the  program  calculates  peak  attenuation  or 
user  defines  peak  attenuations  for  each  target  frequency. 
=  0  original  empirical  attentuation  calculation 

=  1  user  specifies  from  1  to  24  PLAm  values  and  from  1 

to  24  TFm  values 

=  2  least  attenuated  mode  calculations 

=  3  equal  energy  peak  attenuation  calculations 

LGMm  Specifies  whether  program  calculates  peak  attenuation 
using  annular,  rectangular  or  cylindrical  geometry  and 
type  of  input  required. 

=  0  rectangular  or  annular  geometry  with  inputs  EDHm 
and  ELOm  required 

*  1  cylindrical  geometry  with  inputs  RlWm  and  TLm 
required 

(LGMm  =  1  option  only  if  IMAm  =0,  2  or  3) 

The  modified  program  also  keys  on  the  number  of  attenua¬ 
tion  values  input  under  option  IMAm  =1.  If  24  attenuation 


values  are  input,  the  attenuation  at  a  specified  frequency  is 
independent  of  the  attenuation  at  neighboring  frequencies.  If 
less  than  24  attenuation  values  are  input  under  option  IMAm  =  1 
the  program  works  as  originally  designed  where  a  bandwidth  array 
causes  the  attenuation  at  any  frequency  to  be  effected  by  the 
attenuation  at  neighboring  frequencies. 
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SECTION  VII 


PROGRAM  VERIFICATION 


AIRCRAFT  SELECTION 

The  aircraft  selected  for  baseline  noise  studies  and  veri¬ 
fication  of  the  program  modifications  was  the  KC-135B.  The 
standard  version  of  this  aircraft  was  powered  by  four  TF33-P5 
turbofan  engines. 


BASELINE  NOISE  DATA 

With  the  help  of  Lt.  Raymond  Geiselhard,  Mr.  Gregory 
Hughes,  and  Mr.  Robert  Bunting  of  ASD  the  engine  data  required 
to  run  the  noise  prediction  program  for  the  KC-135B  with  four 
TF33-P5  turbojet  engines  was  obtained.  This  data  is  shown  for 
both  takeoff  and  approach  conditions  in  Table  3.  The  noise 
predictions  were  made  with  an  aircraft  Mach  number  of  0.224  at 
takeoff  and  0.140  at  approach. 

Similar  data  was  obtained  for  the  CFM  56  engine.  This 
data,  shown  in  Table  4,  did  not  include  some  of  the  core  noise 
parameters  because  the  core  noise  was  not  thought  to  be  signi¬ 
ficant  for  this  engine. 


COMPARISON  OF  PREDICTIONS  WITH  DATA  AND  CORRECTIONS 

TF33-P5  Noise  Comparison  and  Correction 

The  TF33-P5  engines  on  the  KC-135B  aircraft  was  a  standard 
configuration  tested  for  noise  by  Speakman,  et  al  [16]  with  the 
Aero  Medical  Research  Laboratory  at  WPAFB.  This  data,  shown  in 
Table  5  was  used  to  develop  corrections  to  the  prediction  program 
to  cause  it  to  predict  the  sound  pressure  levels  measured  in 
these  flight  tests. 

The  correction  applied  was  a  function  of  frequency  using 
the  DPB(I)  parameter  from  the  programs  user's  manual. [15]  The 
program  was  corrected  so  that  the  difference  between  the  predic¬ 
tion  and  measured  noise  for  take-off  and  the  same  difference  for 
approach  averaged  to  near  zero.  Actually  this  average  differ¬ 
ence  was  less  than  one  dB  for  20  of  the  24  frequencies.  The 
corrections  were  applied  to  the  jet  noise  portion  of  the  program 
for  frequencies  less  than  1000  Hz  and  to  the  fan  noise  portion 
of  the  program  for  frequencies  of  1000  Hz  or  more.  These 
corrections  are  shown  for  the  TF33-P5  engines  in  Table  5  also. 
The  corrections  are  as  large  as  25.3  dB  near  the  blade  passage 
frequency. 
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TABLE  3 


ENGINE  DATA  FOR  NOISE  PREDICTION 
DATA  FOR  KC-135B  WITH  TF33  P5  ENGINES 

Ambient  Standard  Pressure  59°  F 


1 

Takeoff 

Approach 

Program 

Parameter 

r 

JET  NOISE  DATA  j 

Primary  jet  nozzle  exit  area  -  ft2  i 

3.70 

3.70 

API 

Primary  jet  mass  flow  rate  -  lbm/s 

206.7 

135.6 

WPI 

Primary  exhaust  velocity  returned  j 

to  nozzle  -  ft/s 

1663. 

1020. 

VPI 

Secondary  jet  nozzle  exit  area  -  ft2 

3.65 

3.65 

AS2 

Secondary  jet  mass  flow  rate  -  lba/s  i 

292.1 

235. 

WS2 

Secondary  exhaust  velocity  returned 
to  nozzle  -  ft/s 

1095. 

861. 

VS  2 

Airspeed  -  knots 

250 

160 

FAN  NOISE  DATA 

Pressure  ratio  of  1st  fan  stage 

1.23 

1.064 

FPR45 ( 1 ) 

Pressure  ratio  of  2nd  fan  stage 

1.18 

1.021 

FPR45  (2) 

No.  blades  of  1st  fan  stage 

34 

34 

NB4  5(1) 

No.  blades  of  2nd  fan  stage 

32 

32 

NB4  5(2) 

Fan  rotor  RPM 

6540 

5190 

RNI45 

Rotor/stator  spacing  1st  stage 

53.68 

53.68 

RSS45 (1) 

(a3  percent  of  rotor  blade 
axial  chord) 

Rotor/stator  spacing  2nd  stage 

83.33 

83.33 

RSS45 (2) 

Rotor  diameter  -  ft 

4.184 

4.184 

DIAM4 (1,2) 

Fan  exit  area  -  ft2 

5.29 

5.29 

AREAS (1,2) 

Relative  tio  Mach  No. 

1.17 

.83 

RTS  4  5 

Fan  press,  ratio  at  1.025  tip 

1.6 

1.6 

CFPR4 

CORE  AND  TURBINE  NOISE  DATA 

Combustion  corr.  mass  flow 

rate  -  lbm/s 

30.1 

30.4 

CMF3 

Total  turbine  Dressure  ratio, 

Pi/PQ 

7.04 

6.04 

PP3 

Combustion  exit  temoerature  -  •  R 

2105. 

1720. 

7T3 

Turbine  RPM 

6540. 

5190. 

SS3 

Turbine  outlet  temoerature 

1950. 

1125. 

TU3 

Turbine  tip  dia.  last  stage  -  ft 

3.091 

3.091 

DT3 

No.  blades  in  last  turbine  stage 

80 

80 

BN  3 

Stator/rotor  spacing  -  %  of  stator 
chord 

23.9 

23.9 

CS3 

Primary  mass  flow  -  lbm/s 

208.7 

137.6 

PMF3 
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TABLE  4 

ENGINE  DATA  FOR  NOISE  PREDICTION 
DATA  FOR  KC-135B  WITH  CFM  56  ENGINES 

Ambient  Standard  Pressure  59°  F 


Takeoff 

Approach 

Program 

Parameter 

JET  NOISE  DATA 

Primary  jet  nozzle  exit  area  -  ft2 

4.02 

4.02 

API 

Primary  jet  mass  flow  rate  -  lbm/s 

117.4 

37.7 

WPI 

Primary  exhaust  velocity  returned 

to  nozzle  -  ft/s 

1259.2 

383.2 

VPI 

Secondary  jet  nozzle  exit  area  -  ft2 

3.65 

3.65 

AS2 

Secondary  jet  mass  flow  rate  -  lbm/s 

694.5 

321.0 

WS2 

Secondary  jet  exhaust  velocity 

returned  to  nozzle  -  ft/s 

864.2 

396.6 

VS  2 

Airspeed  -  knots 

250 

160 

FAN  NOISE  DATA 

Pressure  ratio  of  1st  fan  stage 

1.23 

1.064 

FPR4  5(1) 

Pressure  ratio  of  2nd  fan  stage 

1.18 

1.021 

FPR45  (2) 

No.  blades  of  1st  fan  stage 

34 

34 

NB4511) 

No.  blades  of  2nd  fan  stage 

30 

30 

NB45 (2) 

Fan  rotor  RPM 

4594 

2098 

RN145 

Rotor /stator  spacing  1st  stage 

(as  percent  of  rotor  blade  axial 

chord) 

120 

120 

RSS45 (1 ) 

i  Rotor/stator  spacing  2nd  stage 

122 

122 

RSS45 (2) 

Rotor  diameter  -  ft 

4.25,  4.01 

4.25,  4.01 

DIAM4 (1,2) 

Fan  exit  area  -  ft2 

5.35,  5.25 

5.34,  5.25 

AREAS (1,2) 

Relative  tip  Mach  No. 

.913 

.4 

RTS  4  5 

Fan  press,  ratio  at  1.025  tip 

1.3 

1.3 

CFPR4 

Frequency 

(I) 


1000 


3150 

4000 

5000 

6300 

8000 

10000 


TABLE  5 

NOISE  DATA  FOR  THE  KC-135B 
SLANT  RANGE  1000  ft 


Measured 
SPL  Takeoff 


78.3 

73.6 

76.5 

87.5 

91.9 
93.2 

87.1 

90.7 

92.9 

91.1 

90.2 

90.7 

89.2 

88.2 
87.5 
87.2 

87.2 

88.1 

97.2 

90.8 

80.3 
81.0 
74.7 
67.0 


Measured 
SPL  Approach 


69.6 
66.3 

68.7 
75.6 

81.9 

82.8 

76.1 

79.9 

79.2 

80.3 

78.8 

79.5 

77.4 

78.6 

83.1 

78.4 

81.3 

97.5 

91.7 

79.7 

83.9 

77.6 

73.3 

62.1 


Jet  Noise 


Fan  Noise 


Corrections  Corrections 


DPB(I)  db 


14.7 
34.4 

26.7 


-11.5 

-12.7 


-8.5 

-11.1 

-12.5 

-11.0 

-14.0 

-12.8 


DPB(I)  dB 


-10.8 

-25.3 

-14.8 

-7.6 


-10.5 

-20.1 
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CFM  56  Noise  Comparisons  and  Correction 

The  CFM  56  predictions  were  corrected  based  on  static  test 
data  taken  in  a  100  ft  diameter  circle.  This  work  was  done  in 
cooperation  with  Mr.  Ronald  Tagg  of  ASD  as  part  of  a  previous 
Propulsion  Laboratory  contract.  The  correction  was  applied  as  a 
function  of  angle  from  the  engine  centerline.  This  was  done  via 
the  DOB(I)  parameter  of  the  program  users  manual. [15]  The 
resultant  corrections  are  shown  in  Table  6.  The  corrections  are 
as  large  as  11.6  dB  in  the  critical  region  near  90  degrees  from 
the  engine  centerline. 


BASELINE  AND  LINED  NOISE  PREDICTIONS 

The  TF33— P5  engined  KC-135B  data  were  input  to  the  noise 
prediction  program  for  four  duct  lining  conditions.  They  were: 

1.  Hard  walled  ducts. 

2.  Three  foot  lining  optimized  for  the  blade  passage 
frequency  (3488  Hz)  as  inlet  duct. 

3.  Six  foot  lining  optimized  at  3200  Hz  as  inlet  duct. 

4.  Six  foot  lining  with  splitter  optimized  at  3200  Hz  as 
inlet  duct. 

Lining  number  2  above  was  run  using  the  equal  energy 
option.  The  others  were  run  using  the  original  program  duct 
attenuation  predictions.  The  linings  effect  during  approach  was 
much  greater  than  their  effect  during  take-off  because  the 
take-off  was  dominated  by  jet  noise.  The  jet  noise  prediction 
was  not  effected  by  inlet  linings.  The  effects  of  the  linings 
during  approach  is  shown  in  Figure  4.  This  figure  shows  that 
appoach  noise  is  almost  all  at  blade  passage  frequency  for  this 
engine.  The  splitter  lining  caused  about  a  7  dB  noise  reduction. 
The  three  foot  lining  with  no  splitter  had  an  almost  negligible 
efect.  This  is  because  large  ducts  at  this  high  frequency  have 
low  attenuation.  The  splitter,  by  dividing  the  effective  duct 
size  in  half,  causes  the  product  of  duct  size  and  frequency  to  b 
small  enough  that  the  lining  is  more  effective. 

The  CFM  56  engined  KC-135B  data  were  input  to  the  noise 
prediction  program  for  four  duct  lining  conditions  also.  They 
were : 

1.  Hard  walled  ducts. 

2.  Three  foot  linings  optimized  50%  at  2500  Hz  and  50%  at 
3200  Hz  for  both  fan  inlet  and  exhaust. 

3.  Six  foot  linings  optimized  50%  at  1260  Hz  and  50%  at 
2500  Hz  for  both  fan  inlet  and  exhaust. 

4.  Six  foot  linings  optimized  90%  at  1260  Hz  and  10%  at 
2500  Hz  for  both  fan  inlet  and  exhaust. 

Lining  number  2  above  was  run  using  the  equal  energy 
option.  Lining  number  4  was  run  using  the  least  attenuated  mode 
option.  The  others  were  run  using  the  original  program  duct 
attenuation  predictions.  Again  the  effects  of  the  linings  were 
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Figure  4.  TF33-P5  Approach 
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significant  only  during  approach  conditions.  These  effects  are 
shown  in  Figure  5.  This  figure  shows  that  the  fourth  lining  was 
predicted  to  the  the  most  effective  and  that  it  would  be  ex¬ 
pected  to  reduce  the  overall  sound  pressure  level  by  about  8  dB. 
This  lining  was  run  with  the  least  attenuated  mode  option  which 
predicts  less  attenuation  than  the  other  possible  options  at 
these  frequencies.  The  linings  were  more  effective  in  this 
engine  than  for  the  TF33-P5  because  the  dominate  frequencies 
were  lower  which  permits  the  linings  to  have  a  greater  attenua¬ 
tion  . 


PHYSICAL  LINING  CHARACTERISTICS 

The  lining  predictions  for  the  CFM  56  were  based  on  opti¬ 
mized  linings  at  the  fan  blade  passage  frequency  and  its  first 
harmonic.  The  physical  characteristics  of  these  linings  can  be 
obtained  by  utilizing  the  lining  impedance  models  of  Section  II. 
The  perforated  sheet  model  was  used  with  the  assumption  of  a 
seventh  power  boundary  layer  velocity  profile  of  1/2  inch  thick¬ 
ness.  This  thickness  corresponds  to  a  momentum  thickness  of 
about  10%  of  the  total  boundary  layer  thickness.  The  input  to 
the  lining  model  was: 

Perforated  sheet  thickness 
Static  air  temperature 
Static  air  pressure 
Frequency 

Sound  pressure  level 
Flow  Mach  number 

Boundary  layer  momentum  thickness 
Normalized  lining  resistance  R/pc 
Normalized  lining  reactance  X/Pc 

The  resultant  lining  optimized  for  the  blade  passage 
frequency  of  1260  Hz  had  a  depth  of  .218  inches  with  a  percent 
open  area  of  about  5%.  At  2600  Hz  the  optimum  lining  had  a 
depth  of  .192  inches  and  a  percent  open  area  of  about  27%.  A  no 
boundary  layer  analysis  at  this  frequency  gave  a  core  depth  of 
.453  inches  with  about  2%  open  area.  This  sensitivity  of  the 
optimum  lining  open  area  to  boundary  layer  parameters  was  the 
subject  of  Reference  14.  For  these  linings  the  effect  of  the 
linings  on  aircraft  performance  would  be  minimal  because  the 
backing  depths  would  not  significantly  change  the  size  of  the 
engine.  A  procedure  for  computing  these  effects  on  performance 
has  been  developed  in  coordination  with  ASD  personnel.  The  pro¬ 
cedure  was  not  actually  carried  out  for  these  linings  because  of 
the  small  effect  anticipated.  Linings  designed  for  lower  fre¬ 
quency  noise  suppression  would  have  larger  backing  depths  with 
significant  effects  on  engine  size  and  consequently  performance. 


.040  inches 
540  °R 
14.7  psi 
1260  Hz 
160  dB 
-.5 

.05  inches 

5.68 

-8.51 
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SECTION  VIII 

DISCUSSION  AND  RECOMMENDATIONS 


The  lining  impedance  models  of  Section  II  represent  an 
empirical  model  relating  physical  and  acoustic  parameters.  The 
boundary  layer  effects  must  be  considered  as  tentative  until 
theoretical  studies  relating  hole  diameter,  shear  velocity,  and 
other  details  of  the  flow  field  near  the  holes  are  incorporated 
into  the  model.  These  theoretical  studies  have  been  in  progress 
for  years  without  significant  contributions.  The  Section  II 
model  should,  however,  be  adequate  for  a  first  pass  estimate  of 
size  effects  as  required  for  a  good  aircraft  performance  compu¬ 
tation  . 

The  approximate  sound  propagation  parameters  developed  in 
Section  III  permit  a  physically  meaningful  estimate  for  a 
problem  where  a  more  exact  numerical  solution  is  difficult.  It 
is  recommended  that  this  approach  be  extended  to:  (1)  boundary 
layer  effects  via  an  impedance  change  across  the  boundary  layer 
as  approximated  by  an  Eversman  and  Beckemeyer  [17]  type  expan¬ 
sion,  and  (2)  variable  area  ducts  where  ray  acoustics  might  be  a 
natural  way  to  approximate  the  area  variation  effects. 

Sheared  flow  orthogonality  proved  to  be  an  area  where 
considerable  effort  in  seeking  perturbation  solutions  to  the 
governing  equations  produced  no  results.  This  interaction  of 
the  flow  field  with  the  acoustics  is  a  field  of  study  which  is 
related  to  basic  research  in  turbulence  and  represents  an 
important  problem  area.  It  is  recommended  that  numerical 
solutions  of  these  partial  differential  equations  be  used  to 
study  this  problem  further. 

The  study  of  sound  propagation  in  variable  area  ducts 
(Section  V),  produced  some  satisfying  developments.  The  capa¬ 
bility  to  prescribe  any  cylindrical  or  annular  axisymmetric 
geometry  along  with  any  axial  variation  of  wall  impedance  per¬ 
mits  an  impressive  range  of  ducts  to  be  analyzed.  The  bullet 
nosed  inlet  is  a  good  example  of  a  geometry  of  practical  impor¬ 
tance.  It  is  recommended  that  this  numerical  technique  be 
refined  and  used  in  a  parametric  study  involving  axial  impedance 
changes  in  wall  impedance. 

With  flow  included  the  governing  partial  differential 
equations  do  not  permit  an  easy  reduction  to  a  single  variable 
as  in  the  no  flow  variable  area  duct  case.  Difficulty  with 
instabilities  as  time  was  marched  forward  made  progress  diffi¬ 
cult  for  the  numerical  solution  attempted.  Although  stable 
techniques  were  found  they  did  not  permit  tne  range  of  variables 
sought.  It  is  recommended  that:  (1)  Apparent  problems  in 
prescribing  the  time  varying  inlet  condition  be  investigated 
further.  This  apparent  problem  can  even  be  seen  in  the  no  flow 
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results  for  the  bullet  nose  inlet  where  the  first  interior 
imaginary  pressure  seems  to  be  in  error.  For  the  no  flow  case 
the  system  stability  was  enough  to  overcome  this  problem  and 
permit  the  solution  to  approach  a  steady-state  value.  For  the 
flow  case  this  problem  may  be  causing  the  instabilities. 

(2)  Extend  the  flow  solution  which  was  well  behaved  to  include 
compressibility  effects.  This  improvement  is  seen  as  only  a 
minor  change  in  a  portion  of  the  problem  which  is  giving  no 
difficulty. 

The  noise/performance  model  integration  incorporated  the 
basis  axially  uniform  lining  results.  It  is  recommended:  (1)  to 
incorporate  the  results  of  a  multi -segmented  study  as  another 
option,  (2)  to  pursue  studies  of  multi -segmented  variable  area 
ducts  as  a  long  range  goal,  and  (3)  to  continue  the  policy  of  a 
separate  noise  model  with  tools  to  convert  its  output  to  perfor¬ 
mance  related  input.  This  permits  the  flexibility  to  meet 
varying  needs. 

The  program  verification  consists  of  a  verification  that 
the  program  modifications  work  and  of  a  study  of  possible 
acoustic  implications  of  lining  either  the  TF33-P5  or  the 
CFM  56  engined  KC-135.  The  studies  showed  that  the  approach 
noise  could  be  reduced  significantly  with  the  CFM  56  engine  by 
duct  linings.  The  TF33-P5  engine  noise,  however,  was  difficult 
to  attenuate  because  of  its  high  frequency.  With  the  best 
lining  tried  the  CFM  56  engined  aircraft  was  more  than  30  dB 
quieter  than  the  TF33-P5  engined  aircraft  on  approach.  This  is 
with  a  splitter  in  the  TF33-P5  inlet. 
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